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Abstract 

We implement matrix-element corrections to HERWIG parton shower 
simulations for Standard Model Higgs boson production at hadron colliders. 
We study the Higgs transverse momentum distribution and compare different 
versions of HERWIG and resummed calculations. The HERWIG results 
exhibit a remarkable improvement as many more events are generated at 
large transverse momentum after the inclusion of matrix-element corrections. 



1 Introduction 



The Standard Model (SM) of electroweak interactions predicts the existence of the 
Higgs boson, which is responsible for the mechanism of mass generation. However, such 
a particle has not been experimentally discovered yet. Searches for the Higgs boson 
are one of the main goals of the current experiments at the Tcvatron accelerator and, 
ultimately, at future ones, like the Large Hadron Collider (LHC). 

In order to perform such searches, precise QCD calculations are mandatory. Other 
mechanisms may turn out to be useful at hadron colliders [1], yet Higgs production 
via parton fusion is numerically the most important one, especially at the LHC. Here, 
the leading-order processes in the strong coupling constant are qq ^ H and gg — > H, 
with the gluon-gluon fusion mechanism overwhelming the quark-antiquark annihilation 
channel. 

In order to investigate the phenomenology of the Higgs boson, fixed-order calcula- 
tions may be reliable as long as one only considers inclusive observables, such as total 
cross sections; for less inclusive quantities, one needs to account for multi-parton radi- 
ation in order to perform trustworthy phenomenological analyses [2]. Standard Monte 
Carlo (MC) algorithms [3,4] describe parton radiation in the soft or collinear approxima- 
tion, but can have regions of the phase space, so-called 'dead zones', where no radiation 
is allowed. 

In the dead zone, one can however rely on the fixed-order result, as in this region 
the radiation is neither soft nor collinear enhanced. Several methods have recently been 
suggested in order to match parton showers and fixed-order matrix elements [5-8] . In this 
paper we follow the strategy which has been already used to implement matrix-element 
corrections to the HERWIG event generator [3] for several processes: e+e" annihilation 
into quark pairs [9], Deep Inelastic Scattering (DIS) [10], top quark decay [11] and 
vector boson hadroproduction [12]. The dead zone is here populated by using the exact 
next-to-leading order (NLO) tree-level matrix element result and the shower in the 
already-populated region is corrected using the exact amplitude any time an emission 
is capable of being the hardest so far. 

In Section 2 we review the HERWIG parton shower algorithm for the initial-state 
radiation. In Sections 3 and 4 we discuss the implementation of hard and soft matrix- 
element corrections, respectively. In Section 5 we present results on the Higgs transverse 
momentum distributions using different versions of HERWIG and also partonic calcu- 
lations from the literature. In Section 6 we summarize the main results of our work. 
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2 The HERWIG parton shower algorithm 



HERWIG simulates the initial-state radiation in hadron collisions according to a 'back- 
ward evolution', in which the scale is reduced away from the hard vertex and traces the 
hard-scattering partons back into the incoming hadrons [13]. The algorithm relies on 
the universality of the elementary branching probability in the soft or coUinear approx- 
imation. The probability of the emission of an additional soft /collinear parton from a 
parton i is given by: 

The HERWIG ordering variable is = E'^^i, where E is the energy of the parton 
that splits and C,i = with p and pi being the four-momenta of the splitting and 
of the emitted parton, respectively; Zi is the energy fraction of the outgoing space-like 
parton with respect to the incoming one; Pab{z) is the Altarelli-Parisi splitting function 
for a parton a evolving in b. In the massless approximation, = 1 — cos 6', where 6 
is the emission angle to the incoming-hadron direction. For soft emission {Ei -C E), 
ordering according to corresponds to angular ordering. In ((T)) fa{xi, Qf) is the parton 
distribution function for the partons of type a in the initial-state hadron, Xi being the 
parton energy fraction. The quantity 



2ti JqI A;2 Jq^iq^ 



as fQl /•1-Q1/Q2 

dzP[z) 



(2) 



is the Sudakov form factor, expressing the probability of evolution from Q\ to Q\ with 
no resolvable branching. Unitarity dictates that the Sudakov form factor sums up all- 
order virtual and unresolved real contributions. In (^, Qjmax is the maximum value 
of Q, fixed by the hard process, and Qc is the value at which the backward evolution 
is terminated, corresponding, in the case of HERWIG, to a cutoff on the transverse 
momentum of the showering partons. However, since Qc is smaller than the minimum 
scale at which the parton distribution functions are evaluated, an additional cutoff on 
the evolution variable Qi has to be set. If the backward evolution has not resulted in 
a valence quark, an additional non-perturbative parton emission is generated to evolve 
back to a valence quark. Such a valence quark has a Gaussian distribution with respect 
to the non-perturbative intrinsic transverse momentum in the hadron, with a width 
gTint, that is an adjustable parameter and whose default value is zero. 

As the variables and Zi are frame- dependent, one needs to specify the frame where 
the evolution occurs. One can show that, as a result of the Q^-ordering, the maximum Q- 
values of two colour-connected partons i and j having momenta Pi and pj are related via 
= Pi -Pj, which is Lorentz-invariant. For Higgs boson production in hadron 
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collisions, symmetric limits are set in HERWIG: Qimax = Qjmax = ^JVi ' Pj- Furthermore, 
the energy of the parton which initiates the cascade is set to = Qmax = \/Vi " Vy Such 
conditions define the HERWIG frame. It follows that ordering according to implies 
that, in the showering frame, ^ < . 

The region ^ > is therefore a 'dead zone' for the shower evolution. In such a zone 
the physical radiation is not soft or collinear enhanced, but not completely absent as it 
happens in the standard algorithm. It is indeed the purpose of this paper to improve 
the HERWIG parton shower algorithm and populate the dead zone by including matrix- 
element corrections. 



3 Hard matrix-element corrections 

The Born processes leading to Higgs production via parton fusion at hadron colliders are 
gg ^ H (gluon-gluon fusion), of 0(q;|q;), and qq ^ H (quark- ant iquark annihilation), 
of 0{a). In the SM, Higgs production via gluon-gluon fusion is mediated by a quark 
loop. Hereafter, we shall consider only the top quark contribution in the loop, which is 
largely dominant in the SM, with finite top mass. 

Gluon-gluon fusion receives next-to-leading order (NLO) corrections of 0{a'ga) due 
to the elementary processes gg gH, gq qH, qg qH and qq gH. The 
corrections to qq ^ H are instead qq gH, qg qH and gq — * qH^, of 0{asQ:). 

The matrix elements squared for the corrections to gg ^ H can be found in [14, 15], 
where top mass effects are fully included, and in [16], where the authors have considered 
the infinite top mass approximation in the loop. HERWIG uses the formulae of Ref. [14]. 
The actual expressions for such amplitudes are rather involved once the top mass is fully 
taken into account; therefore, we do not report them here for the sake of brevity. The 
real NLO corrections to qq ^ H are instead rather straightforward: the formulae we use 
can be read from Eq. (3.62) of [17], with appropriate Yukawa couplings and crossing. 

In order to implement matrix-element corrections to Higgs production, we follow the 
prescription of Ref. [18], where a method to include higher-order corrections from real 
radiation to a generator of the lowest-order process has been proposed. 

Referring, e.g., to the correction to gluon-gluon fusion given by 

9{Pi)9{P2) ^ 9{p3)H{q), (3) 

^We point out that processes qq gH and qg Hq are NLO corrections to both gluon-gluon 
fusion and quark- antiquark annihilation. However, the actual matrix elements are different whether 
they are corrections to gg — > H, i.e., at one-loop, or to qq H, i.e., at tree-level. 
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we define the partonic Mandelstam variables of process (jH)) 

s = ipi+ P2)^ i = (pi - P3)^ u = {p2- p^f. (4) 

As in [18], we consider a generic Higgs decay H £1^2 and compute tlie differential 
cross sections da^ for tlie 2 — 3 process gg g(-i(-2 and da2 for tlie 2^2 one gg — > ^l^2 
via an intermediate Higgs boson. 

As done in [12], if we assume that the Higgs squared momentum and rapidity in the 
Born process are conserved in the transition from the 2 ^ 2 to the 2-^3 process, we 
find that a factorization formula holds. We obtain: 

where Ai is the matrix element for the process of Eq. Q. In Eq. (0) ctq is the cross 
section of the Born process gg — *• H, fi and /2 are the parton distribution functions of 
the hard-scattering partons in the incoming hadrons, C,i and Xi are the energy fractions 
of the incoming partons in the 2 — 3 and 2-^2 processes, respectively. 

For the gg ^ H processes we evaluate the parton distribution functions and the 
strong coupling constant at m'jj, which is the hard scale of the gg,qq ^ H processes. 
For the gg — > Hg case, we set such scales to the Higgs transverse mass, = <iT + ^"h^ 
qt being the H transverse momentum, which accounts for the additional parton emission 
via matrix-element corrections. Such choices are the same as the ones done for W/Z 
production [12]. 

To include matrix-element corrections, one will have to implement the weight func- 
tion given by Eq. Q. In order to get the phase space which is populated by HERWIG 
and the dead zone, where no radiation is allowed, we can repeat all the steps which have 
been employed for Drell-Yan interactions and report here just the final results obtained 
in [12]. The total phase-space limits, in terms of the variables s and t, read: 

mjj < s < s, (6) 
m]j-s< i < 0, (7) 

where s is the hadronic centre-of-mass energy squared. The value s = mfj corresponds 
to soft radiation, and the lines t = and t = m'jj — s to coUinear emission. 

As shown in Ref. [12], the HERWIG phase space, which corresponds to ^ < z'^ in 
terms of the showering variables, is given by: 

^(7-v^)< s <s (8) 
Lin < i <-'^{3-\ll + ^ 1 , (9) 
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Figure 1: Total (solid) and HERWIG (dashed) phase space for a Higgs mass uih = 115 GeV 
and centre-of-mass energy ^/s = 300 GeV. 



the region which the showering algorithm populates, for niH = 115 GeV, which is the 
HERWIG default value, and ^/s = 300 GeV. The soft and collinear regions are covered 
by the shower and one has an overlapping region, where radiation may come from either 
parton. As expected. Fig. [T] exhibits the presence of a dead zone, where the standard 
HERWIG algorithm allows no radiation. 

Similar arguments are also valid for the other corrections to the gg H Born 
process. Nevertheless, a few comments are in order. The processes qg qH and 
gq — qH are not soft divergent, but only collinear. This implies that the lower limit 
for, e.g., gq qH is tmin = ~ s and there is no overlapping region. Likewise, 
for qg qH, the upper limit is tmax = 0. The process qq gH, when taken as a 
correction to gg H, cannot be interpreted in the parton-shower language as it is 
neither soft nor collinear singular. However, we included this process as well via matrix- 
element corrections: the dead zone will then be the full physical phase space given by 



In order to implement hard matrix-element corrections, we populate the dead zone 
using the probability distribution given by the exact matrix element, i.e., Eq. (jS}, where 
Ai will have to be the correct amplitude squared of the hard-scattering process. More- 
over, in order to maximize the efficiency of the event generation, the fraction of events 
generated according to the different subprocesses gg, qg, gq or qq is proportional to the 
corresponding cross sections in the so-called 'if + jets' process, where the hard process 
is always one of the corrections to gg —>■ H. 



with i 



mm 



^ — s — tmax- In Fig. [T]it is plotted the physical phase space, along with 



Eqs. (EI)-©- 
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Hard corrections to qq ^ H processes are similar to what discussed for the gluon- 
gluon fusion channel. In particular, the phase-space configuration for qq gH is as in 
Fig. P the one for qg — > qH or gq gH is the same as for the analogous corrections to 
gluon-gluon fusion. 



4 Soft matrix-element corrections 



The implementation of soft-matrix-element corrections can be performed using the gen- 
eral method of Ref. [5]. Any time an emission in the HERWIG phase space is capable 
of being the hardest so far, we use the exact matrix element instead of the standard 
HERWIG algorithm. The hardness of a radiation is measured in terms of the transverse 
momentum of the emitted parton with respect to the emitting one. 

The inclusion of the soft correction is performed by multiplying the parton shower 
branching probability by a factor which is the ratio of the HERWIG to the matrix- 
element distribution. It reads: 

(fa d?a - 
dsat dzd^ 

In Eq. fllO|) . J{s,i; z,^) is the Jacobian factor for the transformation (z,^) — > (s, t). As 
the kinematics for Higgs production are the same as for vector boson production, we 
can just use for the Jacobian factor the result reported in Ref. [12]. 

Before closing this section, we would like to point out that matrix-element corrections 
to Higgs production are differently implemented in the PYTHIA event generator [4]. In 
fact, in PYTHIA the parton shower approximation is used in all the physical phase space 
and the exact matrix element corrects only the first emission, rather than the hardest- 
so-far one. Furthermore, the approximation of a top quark of infinite mass is used in 
PYTHIA to define the ratios of the gg Hg and qg — > Hq to gg ^ H cross sections, 
while the latter does use the complete expressions^ . A comparison of phenomenological 
results obtained running HERWIG and PYTHIA is given in [2,19]. 



5 Transverse momentum distribution 

We would now like to present results for the Higgs transverse momentum distribution 
and investigate the impact of matrix-element corrections. In particular, we wish to 
compare HERWIG results without and with such 

^Also, corrections to the qq subprocesses are not available in PYTHIA. 
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50 100 150 200 

qT (GeV) 

Figure 2: Higgs transverse momentum distribution at the Tevatron, according to HERWIG 
with (sohd) and without (dotted) matrix-element corrections. We have set the Higgs mass to 
mn = 115 GeV. 

corrections as well as the improved HERWIG version with the resummed calcula- 
tion of Ref. [20] and the so-called 'Monte Carlo at next-to-leading order' (MC@NLO) 
implementation of [21]. We shall consider Higgs production at the Tevatron and LHC 
and we shall always assume that the intrinsic transverse momentum is gT,int = 0. Unless 
otherwise stated, we shall use the default parton distribution functions of HERWIG, 
but we can anticipate that the relative effect of matrix-element corrections is basically 
the same, regardless of the chosen set of parton distribution functions. 

In Fig. 121 we consider Higgs production at the Tevatron, i.e., pp collisions at ^/s = 
2 TeV, which is the centre-of-mass energy of the current Run 2. In Fig. IHl we plot 
instead the qt distribution at the LHC, i.e. pp collisions at a/s = 14 GeV. We consider 
the HERWIG prediction with (solid histogram) and without (dotted histogram) matrix- 
element corrections. Beyond qr — iTLn/'i the matrix-element corrected version allows 
for many more events. In fact, one can prove that, within the standard algorithm, 
qr is constrained to be gr < fnn- the events at large gr are therefore generated via 
the same exact amplitude. At small qx, the prediction which includes hard and soft 
corrections displays a suppression. By default, after matrix-element corrections, the 
total normalization is still the same and equal to the LO result. It is therefore reasonable 
that the enhancement at large qr implies a reduction of the number of events which 
are generated at small transverse momentum. This result was already found for W/Z 
production [12] and is visible especially at the LHC. 

In Fig. El we plot the improved HERWIG spectrum (solid line) for the LHC, along 
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Figure 3: As in Fig.H but at the LHC. 
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Figure 4: Comparison of matrix-element corrected HERWIG prediction (solid) and + jets' 
(dotted) at the LHC. Here, qq ^ H processes have been turned off. 
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Figure 5: Comparison of matrix-element corrected HERWIG prediction (solid) and the re- 
summed calculation of [50] (dotted). 




Figure 6: Comparison of matrix-element corrected HERWIG (solid) and MC@NLO (dotted) 
predictions. 
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with the result obtained running the H + jets process (dotted hne). In order for such 
a comparison to be rehable, we have turned the qq H hard process off, as 'if + 
jets' does not currently implement the corrections to quark-antiquark annihilation. We 
use the cutoff gT.min = 30 GeV for the 'H + jets' generation. Fig. HI shows that at 
small qr the two predictions are fairly different, but at large transverse momentum 
they agree. This is a reasonable result, since, after matrix-element corrections, large-gr 
events of both processes are generated via the exact amplitude. This is a check that the 
implementation of the hard correction is reliable. 

Next we compare the new HERWIG version with the resummed partonic calculation 
of Ref . [20] , where the authors resummed soft higher-order contributions to the gg ^ H 
process. In fact, the differential distribution da/dqx for the production of a Higgs 
boson of transverse momentum q^ presents terms ~ 1/q^ In"^ (m^/g|,) which get 
arbitrarily large at small g^, i.e., for q^ -C rriH- The leading logarithms (LL) correspond 
to m = n + 1, the next-to- leading logarithms (NLL) to m = n, the next-to- next-to- 
leading logarithms (NNLL) to m = n — 1. 

The authors of Ref. [20] have resummed such enhanced logarithms to NNLL level at 
small qT and matched them to the NNLO result at large qT in order to obtain a reliable 
result over the full qx range. However, for the sake of comparison with HERWIG, which 
includes leading logarithms and only some subleading terms (see, e.g., the discussion in 
Ref. [22] on the comparison of HERWIG and resummation for W/Z production), we use 
the results of [20] in the NLL approximation, matched to the NLO prediction^. 

In order for such a comparison to be trustworthy, we have to do the same assumptions 
as [20]. We use a Higgs mass value uih = 125 GeV and, once we apply matrix-element 
corrections, in Eq. we evaluate the strong coupling constant and the parton distri- 
bution functions at a scale given by the Higgs mass m|^. As in [20], we also use the 
approximation of a top quark with infinite mass in the loop, which is a user-defined op- 
tion of the HERWIG event generator, and choose the MRST2001 leading-order parton 
distribution functions [23]. We finally turn the Born gg- initiated processes off. 

The result of the comparison is shown in Fig. El The normalization and the small-gy 
behaviour are clearly different. In fact, the total cross section is LO in HERWIG and 
NLO in our approximation of the work in [20]. The discrepancy at small transverse 
momenta is instead due to the different logarithmic accuracy which the two considered 
approaches implement. However, the two curves agree at large qT, where the NLO 
calculation dominates. This is another independent check of the reliability of the im- 
plementation of matrix-element corrections. 

■^We point out that the NLO corrections to gg ^ H for the cross section, of order 0{aga), are 
actuaUy LO for the qr distribution and according to the notation of Ref. [20]. For the sake of self- 
consistency, we stick to our conventions and call NLO the fixed-order calculation to which the resum- 
mation is matched. 
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Finally, we would like to compare the results of standard HERWIG after matrix- 
element corrections with the MC@NLO event generator (version 2.2) [21], which imple- 
ments numerically the method of Ref. [7] to simulate Higgs boson production at hadron 
colliders. As discussed in [7], the MC@NLO approach implements both real and virtual 
corrections to the hard-scattering process, in such a way that predicted observables, like 
total cross sections, are correct to NLO accuracy. Moreover, the MC@NLO showers 
turn soft matrix-element corrections off. 

Version 2.2 of the MC@NLO [21] currently includes only the corrections to Higgs 
production in the gluon-gluon fusion channel, hence we shall again have to turn the 
quark-annihilation processes off for the sake of a reliable comparison. In both matrix- 
element corrected HERWIG and MC@NLO generators we set factorization and renor- 
malization scales for the NLO processes like gg Hg equal to the Higgs transverse 
mass. 

In Fig.iniwe show the result of this comparison, which exhibits similar features to the 
one with the resummed calculation. The normalization and the small-gT behaviour of 
the two curves are different, but nevertheless the large transverse momentum predictions 
are in good agreement. 



6 Conclusions 

We have implemented matrix-element corrections to HERWIG parton shower simula- 
tions for direct Higgs production at hadron colliders. We have considered corrections to 
gg ^ H and qq ^ H and have used the exact tree-level NLO matrix element to pop- 
ulate the HERWIG dead zone of the physical phase space and for every hardest-so-far 
emission in the already-populated region. 

We have considered the Higgs transverse momentum qT distribution and have com- 
pared HERWIG predictions with and without matrix-element corrections. We have 
found a remarkable effect of such corrections at large gr, as many more events are gen- 
erated here via the exact amplitude. We have compared the matrix-element corrected 
result with the g-r prediction yielded by the HERWIG 'H + jets' process and have found 
agreement at large transverse momentum. 

We have compared HERWIG provided with hard and soft corrections with a recent 
NLL+NLO soft-resummation calculation and, after making consistent choices for the 
Higgs masses and the scales entering the calculation, we have found very good agreement 
in the large-gr range. 

Moreover, the HERWIG implementation with NLO real corrections has feared rather 
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well against the MC@NLO result, as obtained by using both real and virtual NLO correc- 
tions to the hard partonic process. Besides obvious differences in the total normalization 
(which is LO in standard HERWIG and NLO in the MC@NLO approach) and at small 
qti the large-gT spectra agree well, which is another consistency check of the reliable 
inclusion of matrix-element corrections. 

Between the described implementation and the one available within the MC@NLO 
option, we believe that HERWIG is presently a reliable event generator for Higgs pro- 
duction at hadron colliders both at small and large transverse momentum and that the 
currently-available matrix-element corrections will play an important role to perform 
any analysis on Higgs searches at present and future colliders. In particular, the op- 
tion described here may be the most convenient choice for transverse momentum values 
qr > mu- 

Finally, it should be noted that the HERWIG implementation presented in this paper 
includes both a finite value for the quark mass in the loop of the relevant gluon-gluon 
fusion subprocesses and the corrections to the quark-antiquark annihilation channels, 
thus lending itself to a generalisation of the algorithm to the case of the Minimal Su- 
persymmetric Standard Model, work which is currently in progress. 
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